Exotic nuclei with open heavy flavor mesons 
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We propose stable exotic nuclei bound with D and B mesons with respecting heavy quark symme- 
try. We indicate that an approximate degeneracy of D(B) and D*(B*) mesons plays an important 
role, and discuss the stability of DN and BN bound states. We find the binding energies 1.4 MeV 
and 9.4 MeV for each state in the J p = l/2~ with I — channel. We discuss also possible existence 
of exotic nuclei DNN and BNN. 
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Researches of exotic nuclei have been one of the most 
interesting subjects in nuclear physics. Recently there 
have been much progress in studies of proton- and 
neutron-rich nuclei with large isospin EL and hypernuclei 
and kaonic nuclei with strangeness [2], |3j . Exotic nuclei 
are useful to study various aspects of non-perturbative 
QCD, such as exotic hadrons, nuclear force, dense and 
hot matter, and so forth. They also provide us funda- 
mental information for astrophysics. 

For enlargement of our knowledge of exotic nuclei, the 
variety of multi-flavor is now going to be extended to 
heavier flavors of charm and bottom. So far, several stud- 
ies have been advocated for exotic nuclei with charmed 
baryons [H and charmed mesons Q , in which the interac- 
tion is based on SU(4) flavor symmetry as a straightfor- 
ward extension from strangeness to charm. Recently, in- 
spired by the successful application of the chiral dynamics 
with approximate SU(3) chiral symmetry to strangeness 
sector [3j , the extended version with SU(4) chiral symme- 
try has been applied to charm sector Q. However, the 
dynamics would drastically change in the system with 
heavy quarks, since it realizes not only chiral symmetry 
but also a new symmetry, namely a heavy quark symme- 
try 0). This symmetry has been successfully applied to 
heavy flavor hadrons M, 0, [nl , and other exotic hadrons 

mm 

In this work, we investigate exotic nuclei bound with an 
open heavy flavor meson, D or B meson, with respect- 
ing the heavy quark symmetry. This approach would 
provide us new knowledge of exotic nuclei. Such exotic 
nuclei will be experimentally accessible at future high- 
energy hadron facilities such as J-PARC (Japan Proton 
Accelerator Research Complex) and GSI (Gescllschaft fiir 
Schwerionenforschung) [l3[ • 

One of the remarkable features of the heavy quark sym- 
metry is a degeneracy of pseudoscalar and vector mesons 
as seen in small mass splitting between D and D* mesons, 
and B and B* mesons. Therefore, both pseudoscalar and 
vector mesons are considered as fundamental degrees of 



freedom in the dynamics. The picture for heavy quark is 
completely different from the picture for strange quark, 
in which approximate chiral symmetry is realized. In the 
strangeness sector, only K mesons are important due to 
nearly massless Nambu-Goldstone bosons, and hence K* 
mesons are almost irrelevant. 

The interaction between D (D*) mesons and nucle- 
ons is thus qualitatively different from that between K 
mesons and nucleons. In the charm sector, approximate 
degeneracy of D and D* mesons provides both irDD* 
and nD* D* couplings, and it induces one pion exchange 
potential (OPEP) in the t-channel of the DN-D*N and 
D* N-D* N scatterings. In the strangeness sector, the 
absence of K* mesons leads non-existence of OPEP. In- 
stead, the dominant force is provided by the Weinberg- 
Tomozawa type interaction. Therefore, we consider that 
the DN and D*N interaction is concerned to a long dis- 
tance force with a scale of inverse of w meson mass, while 
the KN interaction is a short distance force with a scale 
of inverse of u> and p meson masses. Such a picture will 
be applied to B and B* mesons with more accuracy. 

In this Letter, we discuss the DN and BN states as 
the smallest system for an open heavy flavor meson and 
nucleons. The s- and u-channcls in the DN and BN scat- 
tering would be suppressed, because the former includes a 
pentaquark with a heavy quark Q, and the latter requires 
a creation of a QQ pair. Consequently, the OPEP in the 
t-channel dominates for this system. In contrast, the DN 
and BN states must be more complicated because of the 
quark annihilation process. For these reasons, we concen- 
trate on the DN and BN states here. As the smallest 
nuclei with D and B mesons, we address possible DNN 
and BNN states with baryon number two. 

The heavy quark symmetry with light quark chiral 
symmetry provides the vertex of 7r mesons and open 
heavy flavor mesons P (D or B) and P* (D* or B*) 



CnHH = gtrH a H b -f„j 5 A% a , 



(1) 



where the multiplet field H of P and P* is defined by 



(2) 
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with the velocity v of the mesons jg]. The conjugate field 
is H a = 70^70, and the index a denotes up and down 
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flavors. The axial current is given by A, 



M = 




V2 



-j-d^M with 



(3) 



where f n = 135 MeV is the pion decay constant. The 
coupling constant \g\ = 0.59 for irPP* is determined from 
the observed decay width T = 96 kcV for D* — > Dir 
[3]. The coupling of nP*P*, which is difficult to access 
from experiments, is automatically determined thanks to 
the heavy quark symmetry. Note that the coupling of 
ttPP does not exist due to the parity conservation. The 
coupling constant g for ttBB* would be different from one 
for ttDD* because of 1/wq corrections with the heavy 
quark mass toq [TBI ]. The recent lattice simulation in 
the heavy quark limit suggests the closed value adopted 
above [16| . This would allow us to use the common value 
for D and B. 

In order to derive the OPEP between and N in 
non-relativistic form, we set v = (1,0) in Eq. ([1]), and 
obtain the irPP* and irP*P* vertcxes 



£ttP*P* 



(e-V)(7f-r) + (£*-V)(7f-f)j , (4) 

(5) 



V2f (f-V)(7f-f), 

JTT 



where the polarization vector of P* are defined by e^' = 
(Tl/V2,±vV2,0) and = (0,0,1), tt is the pion 
field, and the spin-one operator T is defined by Ty X = 



is 



y 'M A,)t ei A) . The nNN 



a) 

P* 2 Q N 
s 1/2 



b) 

P* 4 n N 

u 1/2 



n 



P b 1/2 N 



c) 

P* 2 Q N 
b 1/2 



n 



P* s 1/2 N 



P b 1/2 N 



d) 

P* 4 n N 

u 1/2 



P* D 1/2 N 



FIG. 1: The schematic diagrams of each component in the 
potential V 1 / 2 - in Eq. See the text for details. 



respectively, q is the momentum of the propagating pion, 
and m w is the pion mass (replaced to [i in Eq. 0). A_/v 
and Ap are the cut off parameters for the nucleon and 
heavy meson, respectively. Then, C(r; m) and T(r; m) in 
Eqs. ([7]) and ([8]) arc defined as 



C*(r;m) = 



d 3 p 



1 
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vertex is 



d> 



(2tt) 3 <f 
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where is the nucleon field with spin s, mjv is the nu- 
cleon mass, and g 2 NN /4ir = 13.5 is the coupling con- 
stant. From the vertexes (@}, (J5]), and ©, the OPEPs in 
the PN-P*N and P*N-P*N scatterings are given by 



Vpn- 
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Tp-TN, 

(8) 



VZmNf-rr 3 

x T-aC(r;m^) + STT(r;m^) f P -f N , 



respectively. Here Tp and tjv are isospin operators for 
P w and N. We define the operators S E w = 3(e w -r)(a- 
f)-£^-3 and S T = 3(f-f)(a-f)-f -a. Note that the 
potential ([7]) includes the modified mass scale n 2 = m 2 — 
(nip* — nip) 2 because of the different masses mp. and 
rap for P* and P. In order to estimate the size effect of 
the nucleon and heavy meson, we introduce form factors, 
(A 2 N -ml)/(A 2 N + q 2 ) and (A 2 , - m% )/{A 2 P + q 2 ), in the 
momentum space at each vertex for nNN and ttP^*'P*, 



with 5*12(3;) = 3(<?i • x)((T2 ■ x) — o\ ■ <72, and F{q; m) = 
(A 2 N -m 2 )/{A 2 N + q 2 )x(A 2 p -m 2 )/(A 2 p + q 2 ). 

Now we move to discussion about the pWAT bound 
states. In this analysis, we consider four states, whose 
quantum numbers are classified to J p = l/2~ and 3/2 _ 
with isospin 7 = and 1. The l/2~ states with 7 = 
and 1 are superpositions of three states, 2 5i/2 for PN, 
and 2 Si/2 and 4 Pq/2 for P*N with the standard notation 
2S+1 Lj. The 3/2~ states with 1 = and 1 are superpo- 
sitions of four states, 2 D 3 / 2 for PN, and 4 S 3 / 2 , 4 D 3 / 2 , 
and 2 _Ds/2 for P*N. A similar analysis has been done in 
Ref. [l7| , in which they however do not take into account 
the mixing between different spin states, S\i 2 hi PN 
and P*N, and 4 D 1/2 in P*N. 

With these basis, we explicitly represent the potentials 
© and © as 
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b) PN( 2 S 1/2 )-P*N( 4 D 1;2 ) 
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c) P*N( 2 S 1/2 )-P*N( 2 S 1/2 ) 
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FIG. 2: The plots of each component of the potential V 1 / 2 - with 1 = in Eq. See the text for details 
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for J p = l/2~ and 3/2~, respectively, with the abbrevi- 
ation C m = C(r; m) and T m = T(r; m) for m = m n and 
/i. We confirm that the mixing among the different spin 
and angular momentum states is given in the off-diagonal 
components in the matrices given above. Concerning 
the independnent components in V1/2-, schematic dia- 
grams are shown in Fig.[Q a) PN( 2 S 1/2 )-P*N{ 2 S 1/2 ), b) 
PN( 2 S 1/2 )-P*N( i D 1/2 ), c) PN( 2 S 1/2 )-PN( 2 S 1/2 ), and 
d) P*N{ i D 1/2 )-P*N( 4 D 1/2 ). The kinetic terms are 
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for J p = 1/2 and 3/2 , respectively. Here we define 
A = d 2 /dr 2 + (2/r)d/dr and A 2 = A + 6/r 2 , m pM = 
totv Wp(») /(mjv+mp(»)), and Ampp. ~mp-^mp. The 
eigenvalue equation with the given hamiltonian, iJ pp = 
-ftTjp + Vjp with J p = l/2~ and 3/2 ~, is numerically 
solved by a variational method. The binding energy is 
realized as a difference from the threshold + mp. 

We fix the cut off parameters A^r and Ap, and the 
sign of the coupling constant g. According to Eq. 



the OPEP between two nuclcons is obtained by 
2 



Vnn : 



(fm^) \ [° ; i-° ; 2^+5 12 (f)T m J] r r r 2 .(15) 



We find that Ajv = 940 MeV reproduces a deuteron state 
with the binding energy 2.2 MeV and the relative radius 
3.7 fm. For Ap, we assume the relation, Ap/ An = tn /rp, 
in terms of matter radii of the nucleon rjv and heavy 
meson rp, because the cut off parameter represents the 
inverse of the size of hadrons. In the constituent quark 
model, the sizes rjq and rp are characterized by the fre- 
quencies lon for nucleon and up for heav y m eson, respec- 
tively, in harmonic oscillator potentials [181 ]. In order to 
evaluate ujn and wp, we make use of the charge radii 
which arc relatively well known in experiments. Thus, 
ljn is determined to reproduce the charge radius with 
subtraction of the pion cloud, which is estimated to be 
around 2/3 of the observed value 0.875 fm [T3|, na mely 
0.58 fm, in the analysis in the chiral quark model [l9(. 
Concerning the charge radii of heavy mesons, we use 0.43 
fm for D+ and 0.62 fm for B+ With the constituent 
quark masses m Ut d = 300 MeV, m c = 1500 MeV, and 
ra fc = 4700 MeV, we obtain ^ = 420 MeV and w P = 330 
MeV. Finally, we obtain the ratios t^Itd = 1.35 and 
tn/tb = 1-29, and then Ap, = 1266 MeV and Ap = 1213 
MeV. As for the sign of g, we use g = —0.59. If the op- 
posite sign is used, the binding energy becomes slightly 
larger by a few MeV. From the point of view of G-parity, 
such states would be assigned to DN and BN . 

With these parameters, we obtain the potentials V\/ 2 - 
and V 3 / 2 - in Eqs. (jTTJ) and (|12[) . For example, we plot 
each component in the potential Vi/ 2 - with / = in 
Fig. [2 Each figure corresponds to the diagrams in Fig. [1] 
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TABLE I: The properties of the DN and BiV bound states 
for J p = l/2- with 7 = 0. 
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FIG. 3: The wave functions of the DN and BN bound states 
for J p = l/2~ with 1 = 0. 

The narrow (thick) curves show the D^N (B^*>N) po- 
tentials, and the solid (dashed) curves show the cases 
with (without) the form factor. 

By using these potentials, we find the DN and BN 
bound state solutions for J p = 1/2" with 1 = 0. Their 
binding energies and the relative radii are 1.4 MeV and 
9.4 MeV, and 3.8 fm and 1.7 fm, respectively, which are 
summarized in Table. [I] The wave functions are shown 
in Fig. [3j where the solid, dashed, and dotted curves 
represent PN( 2 S 1/2 ), P*N( 2 S 1/2 ), and P*N( 4 D 1/2 ), re- 
spectively. The narrow (thick) curves show the DN 
(BN) states. We indicate that the P*N( i D 1/2 ) com- 
ponent is crucially important, though it is smaller than 
the PN( 2 Si/ 2 ) as shown in Fig. [3] The binding en- 
ergy is induced mainly by the tensor coupling between 
PN{ 2 S 1/2 ) and P*N( 4 Dj, 2 ), which are in off-diagonal 
terms in V1/2- in Eq. (jlip . Indeed, without this ten- 
sor coupling, we cannot obtain any bound state. For 
more details, one can see the PN( 2 Si/ 2 )-PN( 2 S-l/ 2 ) and 
P*N( i D 1/2 )-P*N( 4 D 1/2 ) components which are diago- 
nal terms in the potential Vi/ 2 - are repulsive as shown 
in Fig. ^p) and d). Obviously, the bound states are not 
realized only by these terms. However, the PN( 2 Si/ 2 )- 
P* N( 4 Di/ 2 ) component with the tensor coupling which 
is an off-diagonal term in V\j 2 - is stronger than the di- 
agonal ones. As a result, the off-diagonal component, 
causing the mixing of the PN( 2 S 1/2 ) and P*N( 4 D 1/2 ) 
states, induces the attraction. 

Thus, the PN-P*N mixing plays an essential role in 
these systems. We see that the BN state is more deeply 
bound and compact than the DN state. This is be- 
cause the the smaller mass splitting between B and B* 
strengthens the BN-B*N mixing, and the kinetic energy 
of B mesons is smaller. In the heavy quark limit, the 
complete mass degeneracy of P and P* induces the ideal 
PN-P*N mixing, which would give the maximal binding 
energy. 

It is worth to emphasize that we can not find DN and 
BN bound states for other channels, J p = l/2~ with 1 = 
1, and J p = 3/2~ with 1 = or 1. As a result, we conclude 



that J p = l/2- with 1 = are the most promising channel 
for detecting stable DN and BN bound states. 

In experiments, the DN and BN bound states would 
be searched in e + e~ collisions, or anti-proton beam with 
deuteron targets [U, ■ Because they do not decay in 
strong interaction, they are experimentally well accessi- 
ble despite of their small binding energies. The DN and 
BN wave functions have the following components, 

\DN) = c° {\D-p)-\D a n))+c? (\D*-p)-\D*°n}) ,(16) 
\BN) = (\B°p)-\B+n))+cf (\B*°p)-\B*+n)) ,(17) 

with some coefficients cf (i = 0, 1 and H = D,B). The 
weak decay processes, D~p — > K + it~it~ +p and B°p — > 
D~n + +p, in the first component would be available for 
reconstruction of the invariant mass. 

We mention that, contrary to the D(B)N states, the 
D(B)N states is more complicated. Although the u- 
channel would be suppressed due to a heavy QQ pair cre- 
ation, the s-channel is not necessarily small. Concerning 
the Born term, it is phenomenologically known in KN 
scatterings that the s-channel plays a minor role [23| . 
However, the A*^ and ir £ c ({,) states below thresholds 

may strongly couple to the D(B)N states. Therefore, 
the coupled channel effects should be included. We note 
that the D(B)N states discussed here are not affected by 
the presence of A*^ nor 7rS c ( b ). 

As further theoretical studies, in addition to OPEP, 
we may include chiral loops, scalar and vector meson 
exchanges, and direct quark exchanges. However, we 
consider that the ranges of their interactions would be 
smaller than the sizes of the DN and BN states, hence 
they would be minor. Nevertheless, the vector meson ex- 
change is interesting because it is compatible with the 
Weinberg- Tomozawa type interaction in chiral symme- 
try [24[ ■ More information about the interactions will be 
provided also from the analysis of molecular pict ure of 
the recently observed exotic charmed mesons [ll|, such 
as X(3872) [H, Z ± (4430) [H], and so forth. 

The stability of DN and BN states should be com- 
pared with the instability of KN states as candidate of 
pentaquarks. It has been theoretically studied that the 
K meson is not bound with a nucleon by the Weinberg- 
Tomozawa interaction with chiral symmetry [27| . In con- 
trast, our analysis shows the D and B mesons are bound 
with a nucleon by OPEP with respecting heavy quark 
symmetry. Therefore, the DN and BN bound states 
may open another way to search for pentaquark with 
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heavy flavors (28j. 

The work on DN bound states is extended to exotic 
nuclei containing D mesons. One can expect that the 
binding energy of the D meson in such exotic nuclei be- 
comes larger as baryon number increases. Such exotic 
nuclei will be investigated in J-PARC and GSI by us- 
ing anti-proton beams with targets of nuclei (2|| [30| . 
Let us discuss the DNN states with baryon num- 
ber two. The possible states for 1—1/2 arc classi- 
fied to \D(NN) 01 ) + \D*(NN_) 10 ) with J p = 0" and 
\D(NN) h0 ) + \D*\NN) 1 .o) + \D*'(NN)o,i) with J p = 1". 
For I = 3/2, they are \D(NN) ,i) with J p = 0" and 
|£)*(7V7V)o j i) with J p = l~, where subscripts denote the 
spin and isospin of nucleon pairs. Only the former two 
states for I = 1/2 contain the attractive DN pairs for 
J p = 1/2- with 1 = 0. Therefore, one can expect that 
the DNN states, as well as the BNN states, for J p = (T 
and/or 1~ with 7=1/2 would be stable. 

It is valuable to compare the DNN and BNN states 
with the K~pp states (K(NN) ,i in our notation) [32l |. 
The K~pp state is considered to have large binding en- 
ergy 20-70 McV. However, the K~pp state decays to 
ttYN by a strong interaction with the decay width 40-70 
McV. On the other hand, the DNN and BNN states are 
stable in the strong decay as discussed above. Therefore, 
the DNN and BNN states may provide more precise 
information about exotic nuclei. More quantitative anal- 



yses including few-body calculation will be presented in 
Ref. |3. 

The charmed and bottom nuclei like DNN and BNN 
states are comparable with the hypothetical existence of 
the kaonic nuclei, in which K mesons are bound in nuclei 
[3l| . The existence of such nuclei may realize if there is 
a sufficient number of nucleons around K meson, though 
KN interaction itself seems not attractive to form KN 
molecule [27| ■ We stress again that in our analysis the D 
and B mesons have sufficiently large attraction to form 
the stable charmed and bottom nuclei beginning with 
baryon number one. 

In summary, bound states of nucleon and an open 
heavy flavor meson are discussed with respecting the 
heavy quark symmetry. It is found DN and BN bound 
states with binding energies 1.4 MeV and 9.4 MeV, re- 
spectively, for J p = l/2~ with 1 = 0, and no bound states 
in other channels. These states are stable in the strong 
decay, and can be observed in the weak decay processes 
DN -> K + n~ir~ +p, and BN -> D~n + + p. The ex- 
istence of DN and BN bound states would provide an 
opportunity to probe new exotic states near the thresh- 
olds, and open a new way to investigate for exotic nuclei 
with variety of multi-flavor explored at future hadron fa- 
cilities such as J-PARC and GSI. 

The authors thank A. Dote for fruitful discussions. 
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